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Abstract 

We report the results of detailed thermodynamic investigations of the Sn2o cluster using density- 
functional molecular dynamics. These simulations have been performed over a temperature range 
of 150 to 3000 K, with a total simulation time of order 1 ns. The prolate ground state and 
low-lying isomers consist of two tricapped trigonal prism (TTP) units stacked end to end. The 
ionic specific heat, calculated via a multihistogram fit, shows a small peak around 500 K and a 
shoulder around 850 K. The main peak occurs around 1200 K, about 700 K higher than the bulk 
melting temperature, but significantly lower than that for Snio- The main peak is accompanied 
by a sharp change in the prolate shape of the cluster due to the fusion of the two TTP units to 
form a compact, near spherical structure with a diffusive liquidlike ionic motion. The small peak 
at 500 K is associated with rearrangement processes within the TTP units, while the shoulder at 
850 K corresponds to distortion of at least one TTP unit, preserving the overall prolate shape of 
the cluster. At all temperatures observed, the bonding remains covalent. 
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I. INTRODUCTION 



Clusters of the semiconducting Group IV elements are interesting not only for their 
potential technological applications, but also from a fundamental point of view. Small 
clusters of Si, Ge, and Sn have an unusual size-dependent structural evolution. Ion mobility 
measurements^! revealed that for small sizes the clusters are prolate, with aspect ratios 
up to three, but undergo a rearrangement to a more compact, spherical form at a size of 
a few tens. The rearrangement is rather abrupt for Si and Ge, but more gradual for Sn. 
Density-functional theory (DFT) studies of ground-state structures&i'Bi'i00 and optical 
response^ confirmed the prolate shape of the small sizes and showed further that many of 
the ground-state structures in this size range contain the tricapped trigonal prism (TTP) 
unit, a nine- atom unit consisting of a triangular prism capped by one atom on each of the 
three side faces. 

Recently, interest was focused on the thermodynamics of the prolate Sn clusters after 
Shvartsburg and Jarroldll! found evidence, based on measurements of ionic mobility, that 
these clusters remain solidlike at temperatures up to at least 50 K higher than the melting 
point of bulk Sn [T m (bulk) = 505 K]. This is the only example so far observed of a small par- 
ticle melting at a higher temperature than the bulk. It contradicts the standard paradigm, 
based on thermodynamic arguments, that a cluster should melt at a lower temperature than 
the bulk because of the effect of the cluster surface .0'IiIElI Many examples of a melting-point 
depression for clusters of mesoscopic size (greater than a few thousand atoms) have been 
demonstrated experimentallyllllll and recently a depression was also observed for small Na 
clusters in the size range of 50 to 200.EH 

There have been a few DFT studies that have implications for the thermodynamics of Sn 
clusters. Lu et a/.0 carried out simulated-annealing studies of Snjy clusters (N < 13), within 
the local-density approximation (LDA), in which the clusters were cooled from a random 
atomic configuration at 1980 K down to K in 22 ps. By following the root-mean-square 
atomic displacements as a function of time, they found that the final ground-state geometries 
were reached roughly half way through the annealing process, at about 1000 K, for 10 < 
N < 13. Majumder et a/.0 calculated energies of some low-lying isomeric structures within 
the generalized gradient approximation (GGA) for the exchange-correlation functional, and 
pointed out that the binding energy per atom was unusually high, only about 11% less than 
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that of the bulk. Finally, we studied the thermodynamics of the Snio cluster the smallest 
Sn cluster whose ground state displays the TTP unit. Using DFT in the LDA, we found the 
TTP unit to break up for temperatures of about 1600 K or higher, yielding a broad peak in 
the canonical specific heat around 2300 K. 

In this paper, we revisit the thermodynamics of Sn clusters in greater detail. Prolate 
structures are identified in ion mobility experiments by having a smaller mobility than a 
spherical structure. Shvartsburg and Jarrold0 observed that the mobility of SnJ- clusters 
(10 < iV < 30) remained small even at 50 K above the bulk melting point. They argued 
that, if the clusters had melted, they would have collapsed to a more compact, spherical 
form, observable by an increase in mobility, contrary to the actual observation. Now, the 
cluster Snio, which we studied earlier,0 contains just one TTP unit and already has a 
compact structure both in its ground state and in high-temperature liquidlike states. Thus 
one would not expect to be able to detect its melting in a mobility experiment. We have 
therefore extended our methods to treat a cluster of twice the size, Sn 2 o- In its ground state, 
Sn2o has a prolate structure, with an aspect ratio of about 2.5, consisting of two TTP units 
stacked end to end. In this paper we study the structural rearrangements, caloric curve, and 
electronic properties of this prolate cluster, and contrast the results with those of Sn 10 . 

As in our earlier paper ,@ we use ab initio Born-Oppenheimer molecular dynamics, within 
the LDA of DFT, and norm-conserving pseudopotentials. This first-principles treatment of 
electronic structure is highly desirable, since the directional covalent bonding in small Sn 
clusters would be hard to describe reliably by a parametrized potential, simultaneously for 
the ground state and for the higher-lying isomers and liquidlike forms crucial for under- 
standing the melting properties. Despite the computational expense of ab initio approaches, 
we obtain a sufficiently long statistical sampling of the ionic phase space as to converge the 
canonical specific heat to of order 10%, as extracted from a multihistogram fit. Our total 
sampling time is of order 1 ns. We believe that this statistical convergence criterion on the 
specific heat gives us a good sampling of the various dynamical and rearrangement processes 
of the cluster function of temperature. 

As we shall see, the Sn 2 o cluster undergoes a number of dynamical processes as a function 
of temperature. One of these is a collapse of the prolate structure toward a more spherical, 
disordered structure. This collapse, which is the point where an increase in ionic mobility 
would be observed, occurs about 700 K above the bulk melting point in our LDA model. 



3 



Thus our calculations confirm, and explain, the experimentally observed stability of the 
prolate form at 50 K above the bulk melting point. El However, at lower temperatures, of 
order the bulk melting point or lower, we also find a large number of dynamical processes 
that lead to a permutational rearrangement of the atoms of the cluster, but which preserve 
the overall prolate shape and would therefore be unobservable in a mobility measurement. 

Since we perform DFT calculations of electronic structure for every ionic configuration in 
our thermal ensembles, we can also study electronic properties as a function of temperature. 
Noting that solid bulk Sn undergoes a phase transition at 286 K between a low-temperature 
covalent form and a higher-temperature metallic form, we examine the bonding character- 
istics of the various isomers of Sn 20 using the electron localization functionlll (ELF). The 
bonding is found to be predominantly covalent at all temperatures considered. However, 
the gap between the highest occupied molecular orbital and the lowest unoccupied molecu- 
lar orbital (HOMO-LUMO gap) is reduced dramatically during the collapse of the prolate 
structure. Finally, although we study only Sn clusters in this paper, we expect that a 
large number of our results would be qualitatively similar for Si and Ge clusters, since the 
structure and bonding of these clusters is similar to those of Sn. 

The paper is organized as follows. In Sec. [TT] we describe briefly our computational and 
statistical approaches. Results for the structural, caloric, and electronic properties are given 
in Sec. [TTI|, and the conclusions are given in Sec. [TV]. 

II. COMPUTATIONAL DETAILS 

All the calculations have been performed using the Kohn-Sham formulation of density- 
functional molecular dynamics (MD),@ using the nonlocal norm-conserving pseudopoten- 
tials of Bachelet et a/.EI in Kleinman-Bylander formll with the Ceperley-Alderll exchange- 
correlation functional. We use a cubic supercell of length 40 a.u. with an energy cutoff of 
16.3 Ry, which is found to provide a convergence of the total electronic energy of better 
than 10~ 6 a.u. To have a sufficiently precise evaluation of the Hellmann-Feynman forces on 
the ions, we ensure that the residual norm of each KS orbital, defined as l-fficsV'i — e ii J i\ 2 
(ej being the eigenvalue corresponding to eigenstate tpi of the KS Hamiltonian i^Ks) ; was 
maintained at 10~ 9 a.u. The ground state and other equilibrium structures have been found 
by carrying out several steepest-descent runs starting from structures chosen periodically 
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from a high-temperature simulation. 

The ionic phase space is sampled by isokinetic MD,!§ in which the kinetic energy is held 
constant using velocity scaling. We split the temperature range 150 K < T < 3000 K into 
about 30 different constant temperatures, performing up to 40 ps of simulation for each. 
We observe no evaporation on the time scale of our simulations, and therefore a reflective 
container (often used with classical potentials) is not explicitly required. Following the 
multiple-histogram (MH) method,0 we then extract the classical ionic density of states Q(E) 
of the system, or equivalently the classical ionic entropy S(E) = kslTiQ(E). With S(E) in 
hand, one can evaluate thermodynamic averages in a variety of ensembles. We focus in this 
work on the ionic specific heat and the caloric curve. In the canonical ensemble, the specific 
heat is defined as usual by C{T) = dU(T)/dT, where U{T) = J Ep(E,T) dE is the average 
total energy, and where the probability of observing an energy E at a temperature T is given 
by the Gibbs distribution p(E,T) = Q(E) exp(—E/kBT)/Z(T), with Z(T) the normalizing 
canonical partition function. In the microcanonical ensemble, the temperature T(E) at 
energy E is defined following the thermodynamic definition as T(E) = [dS(E)/dE]~ 1 . The 
microcanonical specific heat is then defined as C(E) = [dT(E)/dE}~ 1 . More details of the 
extraction of thermodynamic averages can be found in our earlier work.0 

The nature of the bonding can be investigated using the ELF0 along with the charge 
density. Such ELF have been found to be extremely useful for elucidating the bonding 
characteristics of a variety of systems, especially in conjunction with the charge density. For 
a single determinantal wavefunction built from KS orbitals ipi, the ELF is defined as 

X ELF = [l + (D/D h ) 2 ]-\ (1) 

where 

D h = (3/10)(37r 2 ) 5/ y/ 3 , (2) 

D = (1/2) £ |V^| 2 - (V8)| Vp| 2 /P, (3) 

i 

with p = p(r) the valence-electron density. A value of Xelf near 1 represents a perfect 
localization of the valence electron density.0 
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FIG. 1: (a) Some prolate equilibrium geometries having two distinct TTP units. The structure on 
the left is the lowest-energy structure found, (b) Nonprolate equilibrium structures found at high 
temperature T > 1250 K. (c) Structure observed during the transition from prolate to nonprolate 
around T = 1250 K. The excitation energy relative to the ground state is shown as AE, in units 
of eV per atom. 

III. RESULTS AND DISCUSSIONS 

We begin our discussion by presenting some equilibrium geometries of Sn 20 (shown in Fig. 
|l]). We distinguish two broad classes of equilibrium structures: one [Fig. [l|(a)] where all the 
structures have a prolate shape and show the presence of two distinct TTP units, but differ 
in the relative orientation of the TTP units, and a second one [Fig. 0(b)] with nonprolate, 
compact geometries indicating that the TTP units have lost their identity. In Fig. |l](c), we 
also show a typical intermediate structure in which only one of the TTP units has lost its 
identity, but which is still prolate. The most symmetric structure (proposed as the ground 
state in Ref. [II]) is slightly higher in energy than our lowest-energy structure. The only 
difference in these two structures is the relative orientation of the TTP units. Evidently, 
changes in the relative orientation of the TTP units give rise to several low-lying isomers 
[two of which are shown in Fig. 0(a)] having excitation energies of order 0.01 eV per atom 
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FIG. 2: Isosurfaces of the electron localization function xelf for a low-lying prolate form of Sn2o 
containing two TTP units; (a) xelf = 0.7, (b) xelf = 0.55.0 

relative to the ground state. It is also interesting to note that a recent studyii has found the 
most stable structure of Si2o to be very similar to that for Sn2o found in the present work. 

In Fig. 0(b), we have shown two nonprolate equilibrium geometries having a difference in 
binding energy with respect to the ground state of about 0.1 eV per atom. These nonprolate 
structures are somewhat disordered, although some basic units such as tetrahedra, trigonal 
prisms, pentagonal bipyramids, etc., may be observed. The equilibrium structure shown 
in Fig. |I](c) is an intermediate structure observed during the transition from prolate to 
nonprolate. In this geometry, one TTP unit is intact while the other is completely destroyed. 
This structure is 0.075 eV higher in binding energy per atom compared to the ground state. 

It is of considerable interest to investigate the nature of bonding in the Sn 20 cluster. It 
may be recalled that at room temperature bulk Sn exists in a metallic phase (white tin), 
while in our previous workEI we have shown Sn 10 to be a covalently bonded cluster. While 
we therefore expect the bonding within a TTP unit in Sn 2 o to be covalent, it is interesting 
to examine the bonding between the TTP units. In Fig. ^| we show isosurfaces of the ELF 
taken at two different values of Xelf for a typical prolate structure. Recall that a higher 
value of the ELF, approaching unity, indicates a localized distribution of valence electrons. 
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FIG. 3: Canonical specific heat of Sn2o cluster. Co = (3iV — 9/2)&b is the zero-temperature classical 
limit of the rotational plus vibrational canonical specific heatS 

We note that the isosurface for the higher value Xelf = 0.7, shown in Fig. ||(a), connects 
the Sn ions within each TTP unit along the line joining them, implying that the ions within 
each TTP unit are covalently bonded. However, this higher- value isosurface of the ELF does 
not connect the two different TTP units. The bonds between the two TTP units are only 
observed at a rather lower ELF value Xelf ~ 0.55, shown in Fig. 0(b). Evidently the two 
TTP units are weakly bonded to each other compared to the bonding within a TTP unit. 
This leads to a number of low-lying isomers with nearly degenerate energies, of which some 
are shown in Fig. |TJ(a) , and which as noted above differ in the relative orientation of the 
TTP units. 

Turning to the canonical ionic specific heat (shown in Fig. |3]), we note some interesting 
features. First, the main peak in the specific heat, often identified with the "melting" of the 
cluster, occurs around 1200 K, much higher than the bulk melting point (505 K). Second, 
there is also a small peak around 500 K and a shoulder around 800 K. These features can 
be explained by examining the microscopic processes as a function of temperature. At low 
temperatures (below 400 K), the only motion observed is a change in the relative orientation 
of the two TTP units, which spans the class of low-lying isomers [Fig. 11(a)] having two 
distinct TTP units. 

However, at around 500 K, the atoms within each TTP unit start to rearrange themselves 
in such a way that the TTP structure remains intact, but undergoes a reshuffling of at least 
three atoms in every rearrangement, including an interchange of capping atoms with atoms 
of the trigonal prism core. This process is similar to the one observed in Sn 10 S Since this 
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FIG. 4: Normalized potential-energy distributions for 400 K to 1375 K. The distributions are 
obtained from the actual data of an isokinetic sampling run at the kinetic temperature shown, but 
have been smoothed to suppress fluctuations and make the figure clearer. 



rearrangement does not distort the TTP unit, the shape of the cluster remains prolate. The 
distortion of the TTP units themselves begins around 800 K. Starting at this temperature, 
we observe that at least one of the TTP units may be seriously distorted at various points 
in the dynamics, to the extent of losing its identity. The onset of these processes at 500 K 
and 800 K is responsible for the features seen in the canonical specific-heat curve. 

This is also reflected in the potential-energy distributions (PED) for these temperatures. 
In Fig. H, we show normalized PED for a series of temperatures starting at 400 K and 
going up to about 1300 K. Let us recall that the PED p(E, T) at temperature T follows 
the Gibbs distribution p(E,T) = Q(E) exp(— E /ksT)/ Z{T), where Z(T) is the canonical 
partition function and Q(E) is the density of accessible states. An examination of Fig. |] 
reveals a slight relative broadening of the PED around 800 K, indicating that other modes 
of excitation are also available at this temperature, apart from the relative orientation of the 
TTP units and the internal rearrangements within the TTP units that take place at lower 
temperatures. This broadening is reflected in the ionic specific heat as a shoulder around 
800 K. 

As the temperature rises above 800 K, the relative motion between the two TTP units 
becomes more vigorous. The TTP units, apart from undergoing internal distortion, also 
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FIG. 5: The coefficient e pro [see Eq. (Q)] describing the degree of prolate deformation of the cluster 
as a function of temperature. Note the significant decrease around 1250 K. 

start to come sufficiently close to one another as to initiate the exchange of atoms. This 
process ultimately leads at about 1200 K to the collapse of the prolate structure toward a 
more compact structure such as those shown in Fig. |l|(b). In fact, in the temperature range 
1200 to 1300 K, we occasionally observe that the collapsed structure reforms into a prolate 
structure. Thus the time scale at these temperatures for "hopping" between a prolate and 
a more compact structure is of the same order of magnitude as that of our sampling runs, 
namely, about 40 ps. 

To summarize, there are three main modes of ionic motion observed during our simula- 
tions up to 1200 K: (i) changes in the relative orientation of two Sn 10 units, (ii) an internal 
rearrangement of ions within each TTP unit, observed after 500 K, and (hi) the distortion 
of at least one TTP unit accompanied by the interchange of atoms between the two TTP 
units, observed after 800 K, leading ultimately to the fusion of the TTP units and a shape 
change at 1200 K. 

Since it is the collapse of a prolate shape toward a more spherical shape that causes 
a change in the diffusion coefficient in an ion mobility experiment, 111 it is interesting to 
investigate these shape changes more explicitly. In Fig. |5|, we show the thermally averaged 
value of the deformation coefficient e pro , calculated for temperatures ranging from 200 to 
3000 K, and averaged over 35 ps for each temperature. For a given ionic configuration, e pro 
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FIG. 6: The root-mean-square bond-length fluctuation 5 Tms [Eq. (0)] of Sn2o- 



is defined as 
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pro Q 2 + Q 3 ' 

where Q\ > Q 2 > Q3 are the eigenvalues, in descending order, of the quadrupole tensor 

Qij = X! Rij- 
1 



(4) 



(5) 



Here i and j run from 1 to 3, / runs over the number of ions, and Ru is the zth coordinate 
of ion I relative to the cluster center of mass. We note that e pro will be unity for a spherical 
cluster and greater than unity for prolate deformations. At very low temperatures, we see 
that e pro ~ 6, corresponding to the ground-state structure in Fig. |T](a), which can be regarded 
roughly as having a uniaxial prolate deformation with aspect ratio ^e pvo ~ 2.5. This value 
of e pro is maintained until about 800 K, where fluctuations are observed, corresponding to 
the distortion of one or both TTP units, as described above. However, the cluster still 
maintains an overall prolate shape until about 1250 K, when there is a very sharp drop in 
e pro toward a value of about 2 to 3, which is maintained up to temperatures of 3000 K or 
more. This drop corresponds to the fusion of the two TTP units to give the more spherical 
shapes shown in Fig. 0(b). 

In discussions of cluster "melting" phenomena, it is often useful to consider the root- 
mean-square bond-length fluctuation <5 rms , defind as 

x 2 v ((Rjj) t - (R u } 2 t ) 1/2 

0ms iv(iv - 1) e, 



i>j 



(Ruh 



(6) 
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FIG. 7: The canonical (solid line) and micro canonical (dashed) caloric curves for Sn2o- In the 
microcanonical case, the temperature is defined as T(E) = [dS(E)/dE]~ 1 for excitation energy E. 
In the canonical case, the temperature is the heat-bath temperature, and the excitation energy is 
the thermally averaged excitation energy. 

where N is the number of ions in the system, Rjj is the distance between ions / and J, and 
(• • -) t denotes a time average over the entire trajectory. The calculated <5 rms is shown in Fig. ||. 
It can be observed that the traditional Lindemann criterion of 5 Tms = 0.1 signaling bulk 
melting is reached quite early, around 650 K. This happens because <5 rms is sensitive to the 
permutational rearrangement of ions within a TTP unit that occurs at these temperatures. 
This highlights the difficulty of giving a precise definition of "melting" in a small cluster. 
In the temperature range 400 to 1300 K, Sn 2 o undergoes a series of complex rearrangement 
processes. In some sense, an individual TTP unit could be said to have melted already at 
500 K, since after a sufficiently long time the reshuffling motion found at this temperature 
allows ions to permute throughout a TTP unit. After about 800 K, when movement of 
ions between the two TTP units becomes possible, ions can permute throughout the whole 
cluster. However, the dominant peak in the specific heat and the collapse of the overall 
prolate shape of the cluster do not occur until 1200 K. Note that a small step in 5 rms can 
also be observed near 1250 K, corresponding to this collapse. In this final, collapsed state, 
for temperatures of 1300 K and upwards, the ionic motion is diffusive and the cluster has a 
highly liquidlike character. 

The canonical and microcanonical caloric curves are shown in Fig. [7|. The canonical curve 
shows a mild kink extending across the temperature range 800 to 1400 K, corresponding to 
the peak in the canonical specific heat in Fig. |3|. The microcanonical curve has a somewhat 
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FIG. 8: The microcanonical specific heat for Sri2o- For Co, see caption to Fig. ||. 

more distinct step around 1250 K, but the curve does not develop an S bend, as is theoret- 
ically possible in the microcanonical ensemble.0 No S bend was observed either for Sn 10 in 
our previous work using the same LDA model.0 As is well known, S bends for the solidlike 
to liquidlike transition in clusters typically occur only for special sizes, such as the complete 
Mackay icosahedra for the Lennard- Jones potential.0 Corresponding to the absence of an S 
bend in the caloric curve, the microcanonical specific heat C(E), shown in Fig. ||, is every- 
where finite and positive in the energy range studied. The curve C(E) shows the same basic 
features as the canonical curve C(T) in Fig. |3], but with higher resolution. In particular, 
the shoulder at 800 K and the main peak at 1200 K in C(T) become clearly resolved in 
C{E) as a small peak at an excitation energy E = 0.22 eV per atom and a dominant peak 
at 0.36 eV per atom, respectively. This allows an extraction of the effective "latent heat" L 
for the collapse of the prolate form, L w 0.042 eV per atom, taken to be the area under the 
main peak in C(E) and above a line joining the points on the C(E) curve at E = 0.3 and 
E = 0.42 eV per atom. 

Interestingly, the collapse of the prolate form around 1250 K is accompanied by a rather 
drastic reduction in the HOMO-LUMO gap. In Fig. || we show the thermal average (taken 
over 7.5 ps) of the HOMO-LUMO gap as a function of temperature. The gap, which is 
fairly constant at about 0.036 a.u. up to 1200 K, drops sharply by almost a factor of four 
around 1250 K, and continues to decrease as the temperature rises further. At first sight, it 
is tempting to interpret this reduction in the gap as related to the change in the nature of 
bonding from covalent to delocalized metallic. 

To investigate this point in more detail, we have examined the eigenvalue spectrum of a 



13 



0.04 



^0.035 h 
d 

o 

o 0.025 



P 0.02 - 
i— 1 

1 

2 0.015 h 

o 

ffi 0.01 



0.005 



300 600 900 1200 1500 1800 2100 
Temperature (K) 



(a) 



(b) 



FIG. 9: The H0M0-LUM0 gap, averaged over 7.5 ps, as a function of temperature. 
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FIG. 10: The eigenvalue spectrum for (a) the ground state and (b) a nonprolate structure occurring 
at high temperature. The arrow indicates 40th state. 



typical prolate and nonprolate structure, along with an analysis of the number of nearest 
neighbors of these structures. Although the nonprolate structure is more compact, the 
average number of nearest neighbors in both the prolate and nonprolate structures is almost 
the same. However, in the case of nonprolate structures, the bond angles are not optimal, 
resulting in strained bonds. These strained bonds can be considered to introduce disorder, 
which spreads out the energy levels by lifting degeneracies and splitting groupings of energy 



levels. This is evident from the Kohn-Sham eigenvalue spectrum [shown in Fig. |10[ of a 
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FIG. 11: Isosurface of the electron localization function xelf for a randomly chosen high- 
temperature structure of Sn2o- Xelf = 0.7 and the structure is shown in the inset. 

prolate and nonprolate structure. Although the total occupied bandwidth is almost the 
same, the eigenvalue spectrum of the prolate structure shows distinct bands, in contrast 
to the spread-out eigenvalue spectrum of the nonprolate structure, leading to the observed 
substantial reduction in the HOMO-LUMO gap. To demonstrate that this reduction in the 
HOMO-LUMO gap is due to the disorder in the system and not to a change in the nature 
of the bonding, we have also examined the ELF of some randomly chosen high-temperature 
nonprolate structures. An isosurface of the ELF for one of these structures is shown in Fig. 
|TT] at a value Xelf = 0.7. Quite clearly, there is a significant localization of electron density 
between nearest neighbors, and the bonding is covalent. 

Finally, it is interesting to note that, although in Sn2o the dominant peak in the specific 
heat occurs at a much higher temperature (1250 K) than the bulk melting point (505 K), 
the peak for Sn 10 occurs at a still higher temperature (2300 K).0 To understand this, we 
note some features of the two clusters. Snio contains one TTP unit, whereas Sn2o contains 
two. Although the nature of the bonding within each TTP unit is covalent, in Sn 2 o the two 
TTP units are weakly bonded to each other. When ions acquire sufficient kinetic energy, 
they can cross from one TTP unit to the other, and the two TTP units in Sn2o lose their 
separate identity. That is, the presence of another TTP unit provides an extra channel for 
ions to move. In Sn 10 , the main peak in the specific heat is associated with the breaking of 
covalent bonds within the TTP unit. 

In future work, it will be interesting to investigate the transition from a prolate ground 
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state to a nonprolate ground state as a function of cluster size. 



IV. CONCLUSIONS 

The present ab initio simulations on Sn2o, along with earlier work on Snio,^ provide in- 
sights and an explanation of the higher-than-bulk "melting point" observed experimentally.li! 
There are two basic reasons for this behavior: (i) the covalent nature of the Sn-Sn bond in 
the cluster, compared to the metallic bonding of bulk Sn at room temperature, and (ii) the 
special stability of the TTP unit. We also note that the "melting" of a cluster as detected in 
an ion mobility experiment is specifically associated with the collapse of the prolate cluster 
shape, which requires a collapse of the TTP units themselves. Although many rearrange- 
ment processes can occur at lower temperatures, the dominant feature of the specific heat 
(or caloric curve) for both Sn 10 and Sn 2 o is associated with the break up of the TTP struc- 
ture. In the cluster containing one TTP unit, Sn 10 , the TTP unit breaks up at the higher 
temperature, since the only processes that destroy the TTP unit are the ones that break 
the bonds. As the size increases, for example Sn 2 o, containing two TTP units, the collapse 
of the two TTP units leading to a compact cluster is aided by a relative motion between 
them, and occurs at a much lower temperature than the break up of the TTP unit for Sn 10 . 

In a general way, in a prolate cluster there will be a competition between the covalent 
bonding, which tends to favor stacks of TTP units, and the surface energy, which tends to 
favor more spherical structures. Our calculations reveal that for Sn2o the surface energy 
does not become the dominant factor until temperatures of about 1200 K. It is reasonable 
to expect that as the cluster size increases, the temperature for the collapse of the prolate 
structure will reduce further, until eventually the ground state structure itself is no longer 
prolate, as observed experimentally. 
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